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Abstract. Upward through-going muons in the 
Lake Baikal Neutrino Experiment arriving from the 
ecUptic plane have been analyzed using NT200 data 
samples of the years 1998-2002 (1007 live days). We 
derive upper limits on muon fluxes from annihilation 
processes of hypothetical WIMP dark matter parti- 
cles in the center of the Sun. 
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I. Introduction 

Since long, galactic rotation curves indicate the ex- 
istence of non-radiating dark matter. This evidence has 
been hardened also by evidences on larger cosmic scales, 
ending up with the largest scales, surveyed with the 
help of the cosmic microwave background radiation. In 
particular, the data suggest a dominant contribution of 
non-baryonic, non-relativistic (cold) dark matter Among 
the favoured DM candidates are weakly interacting 
massive particles (WIMPs), in particular neutralinos as 
the lightest particles in many minimal supersymmetric 
models. 

One strategy to search for a neutralino signal is to 
search for neutrinos produced in annihilation processes 
inside a large gravitational mass hke the Sun. This 
method is well established since the nineties using 
underground and underwater(ice) neutrino telescopes. 
The Baikal Neutrino Experiment has published results 
of a search for neutrinos from WIMP annihilations in 
the core of the Earth 11], while the analysis with respect 
to the Sun is now presented the first time. We are looking 
for high energy neutrinos from the Sun in excess of the 
expected atmospheric neutrinos. The analysis is based on 



data taken with the NT200 neutrino telescope between 
April 1998 and February 2003. 

II. The Baikal Neutrino Telescope 

Located in Lake Baikal, South-East Siberia, the Neu- 
trino Telescope NT200 is operated underwater at a 
depth of 1.1 km since 1998. The telescope detects 
Cherenkov light from upward and downward-going rel- 
ativistic muons. NT200 consists of 192 optical modules 
(OMs) arranged pair- wise on 8 strings, with 12 pairs per 
string. The height of the detector is 72 m, its diameter is 
42 m. Each OM contains a 37-cm photomultiplier tube 
(PMT). To suppress background from bioluminescence 
and dark noise, the two PMTs of a pair are switched 
in coincidence. Since 2005, the NT200 configuration 
was upgraded by additional 3 strings, each 100 m away 
from the center. This upgraded detector of about 10 
Mton (named NT200H-) serves as a prototype cell for a 
later Gigaton volume detector. The status of the Baikal 
Neutrino Telescope is presented at this conference 12). 

III. Selected data samples 

We have analyzed NT200 data collected between 
April, 1998 and February, 2003, with a total of 1007 
live days. Calibration methods and methods to recon- 
struct muon tracks have been described elsewhere |3], 
[4J, [5|. Our analysis is based on data taken with the 
muon trigger. It requires N^^f > n within 500 ns, where 
hit refers to a pair of OMs coupled in a channel. 
Typically n is set to 3 or 4. The detector response to 
atmospheric neutrinos and muons has been obtained with 
Monte Carlo simulations based on standard codes like 
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The Sun survey with 510 Baikal neutrino events 
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Fig. 1. Distribution of 510 upward muon directions in the ecliptic 
plane, with the Sun at the coordinate origin. 



Fig. 2. Distribution of correlation angles between muons and the 
Sun for data sample A of 2376 events: data (bullets) and measured 
background (histogram). 



CORSIKA iSl and MUM 13, the Bartol atmospheric v 
flux JS) and the neutrino cross-sections from jO). 

The offline filter which requires at least 6 hits on 
at least 3 strings ("6/3") selects about 40% of all 
triggered events. To distinguish upward and downward 
going muons on a one-per-million mis-assignment level, 
a filter with several levels of quality cuts was de- 
veloped for the atmospheric neutrino {Vat) analysis 
ifTOl . The atmospheric muons which have been mis- 
reconstructed as upward-going particles are the main 
source of background in a search for neutrino induced 
upward-going muons. To get the best possible estimator 
for the direction, we use multiple start guesses for the 
minimization. For the final choice of the local minimum 
of we use quality parameters which are not related 
to the time information. At the offline filter level ("6/3") 
the angular resolution (^f - rm.s. mismatch angle) is 
about 14.1° for the i^at-sample. The present analysis 
defines two samples - sample A and sample B - which 
are optimized for the low and high WIMP-mass region, 
respectively. They use further differently tight quality 
cuts, resulting in different background contaminations 
and sligthly different angular resolution. The quality cuts 
are applied to variables like the number of hit channels, 
the probability of fired channels to have been hit or 
not, the actual position of the track with respect to the 
detector centre and /d.o.f.. To improve the signal-to- 
background ratio we used only events with reconstructed 
zenith angle 6 > 100°. This results in 2376 and 510 
upward going muons for sample A and B, respectively 
(with z/Qt-angular resolutions 5' = 5.3° and ^l* — 3.9°). 
Both '3/ values are much bigger than the visible size of 
the Sun. However the angular window for a signal search 
may be even larger, since at least at energies below 100 
GeV the kinematic al angle between neutrino and muon 
dominates. 





-2 


E 


-3 


o 


-4 


_o 


-5 




-6 




-7 




-8 




-9 



100 200 300 400 500 600 700 800 

Ev [GeV] 

Fig. 3. Neutrino effective area of NT200 versus neutrino energy for 
different selection cuts. Curves from top: offline filter level without 
(short dashed) and with correction (long dashed) for the Sun track; 
for event sample B without (dots) and with correction (solid line) for 
the Sun track. 

IV. Sky-plot analysis 

A map of the ecliptic plane centered at the Sun is 
shown in FiglT] with 510 upward going muons. No 
clustering toward the Sun is observed. The distribution 
of the correlation angles between muons and the Sun is 
shown in Fig|2] The dots refer to angles with the real 
position of the Sun, the histogram to angles with "fake 
Suns", defining the expected background behaviour 

No excess is observed, resulting in upper limits on 
number of muons from the Sun. Table I gives the upper 
limits at 90% confidence level (c.l.) for the two selected 
data samples. We give numbers for three values of the 
half cone to the Sun which are used in the dark matter 
analysis. To obtain upper limits on the muon numbers we 
used the recipes for low statistics given in the PDG ifTTl 
and well as the Feldman-Cousins method ifTSl (numbers 
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for FC in parentheses). 

V. Upper limits on a dark matter signal from 
THE Sun 

To derive limits on neutrino and muon fluxes one 
needs the effective area of the detector. Fig[3] shows 
the NT200 area for neutrinos as a function of neutrino 
energy. The expected numbers of muons generated by 
neutrinos is obtained from the neutrino flux of the 
considered source and the effective area A'^fj of the 
neutrino telescope: 

JEth dE^dil 

(1) 

For neutrinos originating from WIMP annihilation, 
the maximum energy Eup is equal to the WIMP mass 
nT'WiMP, with theory suggesting WIMP masses in the 
GeV-TeV range. With increasing ttlwimp more and 
more channels open up - both in annihilation and decay 
processes. The contribution of energetic ("hard") neu- 
trinos can significantly increase when the WIMP mass 
becomes larger than the respective energy threshold for 
the production of heavy secondaries. 

From calculations of the expected number of muons 
Nf^ according to with the DM neutrino spectra 
from ifTSl . we estimated the size of the angular bin 
for which a WIMP signal could be detected with 90% 
probability. The found relation between half open cone 
toward the Sun and WIMP mass is shown in FigH) 
The three angular bins in Table I are those for 90% 
probabihty of signal detection (M>90GeV, M=50GeV 
and M=30GeV, respectively). Further the upper limits 
on muon numbers were scaled to WIMP masses. Table 
I gives the number of events (data, expected background 
and 90% c.l. upper limit). We assumed neutrinos to be 
produced in two annihilation channels, b anti-b (soft 
channel) and W^W^ (hard channel) which define the 
respective maximum energy E^p for neutrinos. The 
resulting upper limits on muon fluxes at 90% c.l. are 
shown in Fig|5] as function of the WIMP mass, for the 
two selected samples and for soft and hard neutrino 
spectra (normalized to Ethr=^ GeV). The upper limit 
on the WIMP-induced neutrino flux from the Sun is 
found to be = 4.46 x 10^"/cm^^ x yr^^ for 
mwiMP — IQQGeV/c^ and for the annihilation channel 
W-W+ (sample B). 

VI. Discussion 

Up to now all operating neutrino telescopes (see 
the talks of the IceCube and ANTARES collaborations 
at this conference) and the past generation (Baksan 
111,1111 MACRO Ha, Super-Kamiokande 1161) re- 
ported no excess of upward going muons from the 
direction of the Sun when compared to the expectation 
for atmospheric neutrinos. This constrains both the an- 
nihilation rate in the Sun and (in concert with direct 
WIMP searches) the WIMP interaction cross section. 
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Fig. 4. Half cones for muon direction toward the Sun where the 
signal from a WIMP of a given mass could be detected with 90% 
probability. 
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Fig. 5. The Lake Baikal NT200 upper limits at 90% c.l. on the muon 
flux from the Sun for a given WIMP mass (sample A - dashed hne, 
sample B - solid line); "hard" neutrino fluxes are the lower two lines, 
"soff the upper two lines. 

Upper limits on muon fluxes from the Sun obtained 
by neutrino telescopes are shown in Fig|6] (adapted from 
lITTl ; aU limits are for 1 GeV threshold). The Baikal 
curve corresponds to hard neutrino spectra and is for 
mwiMP < 370GeV/ (? obtained from sample B, above 
this value from sample A (extrapolated up to 3 TeV). The 
AMANDA-II ini and IceCube ^ limits are shown for 
hard neutrino spectra. 

The presented results are preliminary, and allow to 
estimate the NT200 sensitivity for high energy neutrinos 
from DM annihilation processes in the Sun. 
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TABLE I 

The 90% c.L. upper limits on the number of moon events from the Sun for analysis samples A and B, calculated 

ACCORDING TO PDG AND (IN BRACKETS) FC, RESPECTIVELY. 
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